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I.  INTRODUCTION 


In  this  report  we  consider  the  Initial  value  problem 


u(0,x,y)  -  u0(x,y) 


Here  u,  F,  G,  H  are  m  X  1  vectors,  and 


F(u)  -  f1 (u)  ♦  f2  (u,|£) 


3u, 


G(u)  -  g1 (a)  ♦  82(u,3x) 


H  -  H(u,x,y). 


(1) 

(2) 


(3a) 

(3b) 


(3c) 


We  shall  assume  that  uQ  belongs  to  the  class  of  functions,  D,  which  are  suf¬ 
ficiently  smooth  that  Equations  (1)  and  (2)  have  a  unique,  strong  solution 
u(x,y,t),  which  for  0  i  t  ST  is  in  the  class  Cp  +  1  of  functions  possessing 
continuous  partial  derivatives  D  au  of  order  through  p+1 ,  for  some  p  £  2. 
Thus,  when  Equation  ; '  !  hyperbolic,  u0,T  must  be  restricted  such  that  no 
shock  formation  occjrs  in  C  £  t  £  T. 

1 .  EVOLUTION  OPERATOR 

Unde'’  the  above  assumption,  there  exists  an  operator1  E^'Etr.tT  with 
the  property  that 


’Jr‘*  -  ETUn  ,  0  S  t  S  t  ♦  t  S  T  («0 

where  L'n  -  u!x,y,t),  Un  ■  j(x,y,t*i  i .  Although  the  applications  normally 
call  for  discrete  values  or.  a  space  lattice,  for  convenience  cf  analysis  we 
prefer  x,  y  in  Equation  ! u )  to  be  variable. 


2.  APPROXIMATE  FACTORIZATION 

A  xa Jo~  problem  of  mode-'r.  numerical  analysis  is  the  ::scove-y  :T  cperat c- 

products 


1 


which  to  p-th  order  accuracy  approximate  the  operator  of  Equation  (^);  i.e., 


.n-M 


ETUn  - 


LTUn 


0  (tP+1) 


(5b) 


Here,  we  shall  be  concerned  with  the  case  p  -  2. 

3.  REGIMEN 

For  economy  of  machine  implementation,  the  operators  are  to  be  the 

simplest  possible.  In  practice,  they  provide  second-order  accurate  approxi¬ 
mations  to  solutions  of  certain  equations  which  are  associated  with  Equation 
(1),  through  the  natural  splitting 


av 

3F(v) 

l  ' 

at  ' 

dx 

V  Ja  , 

3w 

■>  m 

3G(w) 

(6b) 

at 

3y 

39 

at 

H(Q,x,y) 

(6c) 

Tne  corresponding  approximation  operators  for  Equation  (6)  we  denote  by 


v 


,  t  n 
L  v 
x 


n+1 

w 


|  *  0(t  * ) 


(7a) 

( 7t ' 


(7c' 


Tne  method  of  operator  splitting  was  originated  by  Peaceman  and 
Rachford,’  in  deriving  a  variant  of  the  alternating  direction  (ADI)  met  or: 
which  ienas  itself  to  tne  use  of  cyclic  parameters  for  accelerating  conver¬ 
gence.  In  seeding  solutions  of  (1),  the  curse  of  dimensionality  ma.,  be 
avoided  througn  splittings  sjch  as  provided  by  Equations  (5_7),  often  witn 
improved  time  step  rest" i ct ions .  Moreover,  advantage  can  be  taken  of  the  long 
and  successful  ni  story  Cf  research  results  concerning  efficient  rjr.eric-.l 
schemes  for  solving  equations  such  as  Equation  (6).  Particularly  of  note  a-- 
the  advantages  of  MacTcrmaer.’ s  method  for  the  Navier-stokes  Equations,’  an: 


certain  higher  order  shock  capturing  schemes  for  the  Euler  Equations. ^  * 

The  purpose  of  this  work  is  to  provide  rigorous  proof,  in  the  general 
nonlinear  case,  of  the  second-order  accuracy  of  a  splitting  considered  by  Mac- 
Cormack.’  There,  MacCormack  justifies  second-order  accuracy  by  means  of  a 
frozen  Jacobian  analysis  and  a  gain  matrix  approach.  Thus,  his  method  rigor¬ 
ously  establishes  the  result  only  in  the  case  of  a  linear  system.  In  addi¬ 
tion,  we  consider  the  problem  of  obtaining  second-order  splittings  for  systems 
characterized  by  presence  of  deri vati ve-free  source  terms,  as  in  axis- 
symmetric  geometries.  Some  discussion  of  the  optimality  of  the  splitting 
approach  is  given. 


II.  THREE-FACTOR,  SECOND-ORDER  ACCURATE  SPLITTINGS 

Strang*  proves  a  result  on  operator  splitting,  which  is  somewhat  more 
general,  but  whose  content  is  essentially  the  following: 

1 .  SPLITTING  THEOREM  1 

Suppose  operato-s  L1,  LT  are  known,  which  provide,  as  in  Equations 

x  y 

(7a, b),  second-order  accurate  updates  for  solutions  of  Equations  (6a, b). 
Then,  either  of  the  composition  operators  defined  by 


Jn+1  -  (L2  Lt  L2)Un 
x  y  x 


(8a) 


and 


,.n- 


(L2  'T 

y 


L  2)Un 
x  v 


:8b. 


p-ovides  a  seconc-o-de>'  accurate,  three-facto^  splitting  for  the  equation 


9u  9F' u /  +  SO1 u , 
3t  *  3x  'ey 


K')', 


COMMENT 


Using  the  met  he  ds  cf  Strang, *  we  can  show  that  no  two-factor  splitting 
which  emplcys  individually  accurate  operators  can,  ove>'  one  step,  yield  a 
second- order  splitting  fc*~  Equation  (9..  Thus,  among  the  class  cf  ope*' at  o^s 
wnicr.  are  second- o-cer  accurate  over  one  step.  Equations  :8a, O',  are  optimal, 
ir,  terms  of  tr.e  numte-  cr  operators  applied. 


.  ^  T  ^  .  t.  ^  Zr  C  t  r  c ’t  cl"*: 


pi  C  C 


s;' 


3 


S3 


equations  such  as 


3u 

at 


C(t ,x ,y ,u,Dau) 


(10) 


where  Dau  are  derivatives  of  arbitrary  order.  The  equation 


C  ■  a  +  b 


(ID 


signifies  an  arbitrary  splitting,  subject  only  to  the  restriction  that  there 
exist  operators  LT,  which  provide  second-order  updates  for  the  equations 

3  D 


(12a) 


and 


3w 

at 


b 


(12b) 


Sore  complications  emerge  when  C  is  an  explicit  function  of  t;  however,  these 
do  not  concern  us,  as  we  shall  not  require  explicit  time  dependence  of  C. 


III.  FOUh-FACTOR  SECOND-ORDER  ACCURATE  SPLITTING 

New,  observe  that  if  the  operator  sequence  of  Equations  (8a, b)  is  applied 
twice,  six  operator  applications  are  necessary  in  order  to  advance  a  2t  tire 
increment.  KacCo^mac'x3  seems  to  have  been  first  to  notice  that  a  more  eco¬ 
nomical  second-order  update,  over  time  increment  2t,  can  be  obtained.  He  con¬ 
siders  cyclical  applications  of  the  operator  sequence 

,,n*  1  , .  t  .  t  ,  t  i  t  \i,n 

U  =(LLLL)U  (Ip, 

x  y  y  x 

His  j ustif i cation  cf  second-order  accuracy  is  sketched  below: 

Consider  a  Fcu^ie,'  mode 


,,,  i(Xx  +  ry) 

U  ( x ,  y ,  t )  -  A  ( t )  e 


(it. 


Ey  applying  the  operator  sequence  in  Equation  (13)  to  the  Fourier  mode  of  Equ¬ 
ation  (14),  with  frozen  Jaccoian  matrices  Jp,  Jq,  it  emerges3  that  Equation 
(13'  produces  a  gain  matrix 


4 


(15) 


gxgygygx 


which  differs  only  by  third  order  terms  from  the  exact  gain  matrix  obtained 
when  Equation  (14)  is  substituted  in  Equation  (9).  Thus,  Equation  (13)  has 
been  rigorously  justified  second-order  accurate  only  for  linear  systems  Equ¬ 
ation  (9),  with  otherwise  locally-linearized  second-order  accuracy. 

However,  by  approximately  factoring  the  full-step  (middle)  operators  in 
Equations  (8a, b),  we  now  show  that  MacCormack's  cyclically  reversed  sequence 
Equation  (13)  is,  in  general,  second-order  accurate,  subject  only  to  the  re¬ 
strictions  required  for  proving  Strang's  splitting  theorem. 

Suppose  Lt  is  second-order  accurate  for  the  equation 


||  “  g(x,y,Daz) 


(16) 


Then,  tc  within  terms  of  third  order,  it  is  required  that 


Ltz  =  z  +  t  g  (x,y,Daz)  +  Ba  •  Dag 


(17) 


where  B  are  the  Jacobian  matrices  of  g  with  respect  to  the  derivatives  D  z. 
Thus,  3 


,  1  ,  .  -t  s 

L  V  L  Z  .1  = 


t  g  (x,y,D  z  ♦  tD  g)  + 


2  v- 

2  a  a 


D°g 


7 

0(f) 


16 


Expanding  the  second  term  ir.  the  right  member,  we  see  that 


Equation  (9). 


IV.  SPLITTING  IN  THE  PRESENCE  OF  SOURCE  TERMS 

In  recent  research  concerning  second-order  accurate  shock-capturing  algo¬ 
rithms  for  the  Euler  equations,  interest  is  focused  upon  the  problem  of  split¬ 
ting  Equation  (1)  for  the  case  in  which  nonzero  H(u,x,y)  is  present  in  Equ¬ 
ation  (1).  Carofano,7  following  MacCormack* s’  results  for  the  two-dimensional 
case,  intuitively  employs  the  splitting 


Un+?  =  (L  T  L  \  TL  TL  TL  T)Un 
x  y  s  s  y  x 


where  L  T  is  second-order  operator  for  Equation  (6c).  We  now  discuss  opti¬ 
mality  For  Equation  (21),  and  rigorously  establish  second-order  accuracy. 

In  view  of  comments  1  and  2,  it  is  unlikely  that  a  three-factor  product 
of  individually  second-order  operators  can  be  found,  which  over  one  step  with 
time  increment,  t,  provides  a  second-order  update  for  Equation  (1).  What  can 
be  done,  rigorously,  is  to  consider  splittings  which  pair  up  any  two  of  the 
quantities  F,  G,  H  against  the  other.  Typically,  the  splitting 


both  provide  f i ve-f actc" ,  best  possible  in  number,  second-order  accurate 
splittings,  over  one  step  of  increment  i ,  for  Equation  (1)  with  source  tern- 
present.  Among  other  possibilities  similarly  obtained,  the  factorization 


can  be  used  to  establish  second-order  accuracy  for  the  Carofano  splitting  of 

Equation  (21),  when  is  replaced  by  t.  Hence,  Equations  (21,  23 -24)  provide 

equivalent  second-order  accurate  splittings  of  Equation  (1)  in  the  presence  of 
source  terms.  Equation  (24),  applied  over  time  increment  2t,  should  be  most 
efficient,  but  at  the  expense  of  cyclically  modifying  the  time  step. 

Our  final  comment  is  that,  in  terms  of  the  optimal  number  of  operators, 
Strang’s  splitting  of  Equation  (8),  over  time  increment  2t,  is  still  one  oper¬ 
ator  evaluation  more  efficient  than  is  the  MacCormack  version  of  Equation 
(13).  However,  in  many  cases,  stability  restrictions  or  special  problem  idio- 
synoracies  may  mandate  other  priorities. 
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ATTN:  Mr.  S.  Walton 


16 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  coaaents/answers  to  the  iteas/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ Date  of  Report _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Coaaent  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


i 

- — — - — - - - - —  i 

4.  How  specifically,  is  the  report  being  used?  (Inforaation  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


S.  Has  the  inforaation  in  this  report  led  to  any  quantitative  savings  as  far 
as  nan-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ 


6.  General  Coanents.  What  do  you  think  should  be  changed  to  iaprove  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


Naae 


CURRENT 

ADDRESS 


Organization 

Address 


City,  State,  Zip 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Naae 


Organization 

ADDRESS 

Address 


City,  State,  Zip 


(Reaove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 
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